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' We explore the large angular scale temperature anisotropics in the cosmic mi- 

cd ■ crowave background (CMB) due to homogeneous local dust-filled voids in a flat 

^ ■ Friedmann-Robertson-Walker universe with a cosmological constant. In comparison 

O ' with the equivalent dust-filled void model in the Einstein-de Sitter background, we 

^ ■ find that the anisotropy for compensated asymptotically expanding local voids can be 

CN ■ larger because second-order effects enhance the linear integrated Sachs-Wolfe (ISW) 

^ ■ effect. However, for local voids that expand sufficiently faster than the asymptotic 

^ ■ velocity of the wall, the second-order effect can suppress the fluctuation due to the 

Csl ■ linear ISW effect. A pair of quasi-linear compensated asymptotic local voids with ra- 

^ ■ dius (2-3) X 10^ Mpc and a matter density contrast 5^ ~ -0.3 can be observed as 

O ■ cold spots with a temperature anisotropy Ar/T ~ 0(10"^) that might help explain the 

^ ■ observed large-angle CMB anomalies. We predict that the associated anisotropy in the 

o ■ 
^: 

^ ■ Subject headings: cosmic microwave background - cosmology - large scale structure 

H ' 1. Introductio 

Recently, there has been mounting evidence that the statistical isotropy in the large-angle 
cosmic microwave background (CMB) anisotropy may be broken (Tegmark et al. 2003; Copi et al. 
2004; Chiang et al. 2004; de Oliveira-Costa et al. 2004; Eriksen et al. 2004; Hansen et al. 2004; 
Larson & Wandelt 2004; Park 2004; Schwarz et al. 2004; Vielva et al. 2004; Cruz et al. 2005, 
2006 ). Although the significance of any one of the findings is at most 3a, the accumulation of 
anomalies hints at the possibility of the need for new physics to be added to the standard scenario 
(e.g., Luminet et al. 2003; Jaffe et al. 2005). 



local Hubble constant in the direction of the voids could be as large as a few percent. 
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To explain the origin of the anomaUes, it has been suggested that the large-angle CMB 
anisotropy is affected by local inhomogeneities (Moffat 2005; Tomita 2005a,2005b; Vale 2005; 
Cooray & Seto 2005; Rakic et al. 2006). However, none of these explanations has succeeded 
in explaining the specific features of the anomalies, namely, the octopole planarity, the alignment 
between the quadrupole (/ = 2) and the octopole (Z = 3) (Tegmark et al. 2003), and the alignment 
between the low-/ multipoles (/ = 2-1-3) with the equinox and the ecliptic plane (Copi et al. 2004). 
For instance, if one applies a model in which the Local Group is falling into the center of the Shap- 
ley supercluster, the discrepancy between the model prediction and the observed data becomes 
even worse (Rakic et al. 2006). 

Inoue & Silk (2006) explored the possibility that the CMB is affected by a small number of 
compensated local dust-filled voids. In the Einstein-de Sitter (EdS) background, we found that 
a pair of voids with radius 3 x 10^ Mpc and density contrast 5 = -0.3 can account for the 
planarity/alignment features in multipoles with / = 2 and 1 = 3. The cold spot in the Galactic 
southern hemisphere, anomalous at the ~ 3(7 level, can also be explained by such a large void at 
z ~ 1. 

In this paper, we study the signature of compensated local dust-filled voids on the CMB in a 
flat Friedmann-Robertson- Walker (FRW) model with a cosmological constant A. The signature of 
voids in the FRW universe has been extensively studied in the literature (Thompson & Vishniac 
1987; Martinez-Gonzalez & Sanz 1990; Panek 1992; Arnau et al. 1993; Meszaros 1996; FuUana 
et al. 1996; Vadas 1998; Griffiths et al. 2003). However, none of the above discussions have 
considered the effect of the cosmological constant which accelerates the cosmic expansion at late 
times. In contrast to the void model in the EdS background, we expect an extra contribution from 
the linear integrated Sachs-Wolfe (ISW) effect due to potential decay at recent epochs z < 1. It 
is of considerable importance to study whether the second-order effect observed in the equivalent 
void model in the EdS background enhances the linear ISW effect or not (Tomita 2005a, 2005b). 
In a similar manner to our previous work (Inoue & Silk 2006), we use the thin-shell approximation 
for describing the wall. In what follows, we assume that we are outside the voids. In section 2, 
we derive analytic formulae for the temperature anisotropy due to homogeneous dust-filled voids 
in the FRW model with a cosmological constant A, and we study the signature of such voids on 
the CMB. In section 3, we explore a model of local dust- filled voids that agrees with the observed 
anomalies on large angular scales. In section 4, we estimate the mean contribution from the local 
dust-filled voids in the cold dark matter (CDM) cosmology. In section 5, we summarize our results 
and discuss some unresolved issues. 
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2. Dust-filled void model 

2.1. Cosmic expansion 

We consider a flat Friedmann-Robertson- Walker (FRW) model with a cosmological constant 
A as the background universe. In what follows, we use units where the light velocity c is nor- 
malized to 1. In spherical coordinates {r^O^cj)), the flat background FRW metric can be written 
as 

ds^ = -dt^+R\t)(dr^ + r^dn^), (1) 

where t is the time, R denotes the scale factor, and dO,^ is the metric for a unit sphere. Next, we 
consider a homogeneous spherical dust-filled void with a density contrast Sm < 0. We assume that 
the size of the void is sufficiently smaller than the Hubble radius H~^. Then, the metric of the 
homogeneous void centered at the origin can be approximately described by the hyperbolic FRW 
metric as 

ds^ = -dt'^+R%')(dr'^+Rlsmh\r' /Rc)dn^), (2) 

where t' is the time, R' denotes the scale factor, and R^ is the comoving curvature radius. The scale 
factor R for the flat FRW background at time ti can be expanded up to third order in At -t2- h in 
terms of the deceleration parameter q2 and the jerk parameter j2 



RR 
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as 
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, where a dot means the time derivative d/dt, H2 = H{t2), and R{t2) = 1. For a FRW universe 
with matter and A, we have q = (Qm-2QA)/2 where Q„ and Qa are the density parameters for the 
matter and the cosmological constant A, and the jerk is j = -(ilm + ^a)- The density parameter and 
the Hubble parameter (denoted by primes) for the void are written in terms of the matter density 
Pm, the matter density contrast 5^ = Spm/pm-, and the Hubble parameter contrast Sh as 



The deceleration and jerk parameters for the void are 



H' = (1 + Sh)H. 



(5) 
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Equations (I3I4I5I6I) give the scale factor R'{t') for the hyperbolic FRW void. From the Friedmann 
equation, the cosmological time for the flat universe that consists of matter and the cosmological 
constant A is 

2 



t = sH- 



s = 



■An 



Vl-^]^+l 



(7) 



2.2. Time solution 

In what follows, we calculate the evolution of the internal time t' for an expanding void with 
a peculiar velocity in terms of the external time t up to order 0[(rv/H~^)^], where the comoving 
radius of the void is expressed as r,,(0 in the external coordinates and as r[,{t') in the internal 
coordinates. To connect the two metrics at the shell, we require the boundary conditions: 

^2 .z,/..^„/..^^2^ .RVfKit'f 



(R(tM)Y = (R'(t')r[.{t')y\^l+ j , (8) 

and 

-dt^+R^{t)dr^ = -dt'^+R'\t')dr'^. (9) 

where we have assumed that / <^ Rc. Up to order 0[(ry///"^)^], the curvature term in equation ([8]) 
is negligible. Then equation ([8]) and Q yield 

dt' = dt (^\ -RRr,r,6H + ^R^r^Sji^ . (10) 

Let us assume that the expansion of the thin shell (=wall) in the external coordinates is expressed 
as ry{t) oc t'^, where /3 is a constant. Equation (flOl ) can be explicitly written as 

dt'=(i-^e+\sj,e)dt, (11) 



^ = ^^^. (12) 



where we define 

H- 

Note that Sh and ^ are functions of t rather than t'. As we shall see in the subsequent analysis, 
the anisotropy turns out to have a leading order of Therefore, omitting the curvature term in 
equation ^ in the process of deriving dt' can be justifiec|l]- Using equations (I3I4I5I6I7I) . and (fT2l) . 
equation (fTT)) can be integrated to yield the finite time difference A?' = t2-t[, 

A^' = ClA^+c2(AO^ + c3(AO^ (13) 

where c, for /= 1 , 2, 3 are functions of {3 and 5// at ? = ?2- 



'The curvature correction yields additional terms of order in equation ( ITTI ). 



-5- 



2.3. Temperature anisotropy 



We denote quantities at the time the photon enters the void and those at the time the photon 
leaves by the subscripts "1" and "2", respectively. Primes denote quantities measured by a co- 
moving observer in the interior coordinate system, and the unprimed quantities are measured by 
a comoving observer in the background universe just out of the shell of the void (see figure 1 in 
Inoue & Silk 2006). 

To calculate the energy loss, we apply two local Lorentz transformations at each void bound- 
ary. The first is to convert the photon four-vector momentum in the comoving frame in the back- 
ground universe to the frame in which the shell is at rest. The second is to convert it to the frame 
in the comoving frame inside the void. 

The four- vector momentum of the photon that enters the void is 



ki=£i 



/ 1 \ 




(14) 



where V'l is the angle between the normal vector of the void shell and the spatial three-vector of 
the momentum of the photon that enters the void. After the photon passed the shell, the four- vector 
is converted to 



= E, 
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\ 



1 \ 
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/ 
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(15) 



(16) 



where vi and v\ are the velocities of the void shell at the time t = t\ and 7-factors are defined as 
7i = 1/(1 -Vj)^/^ and 7^ = 1/(1 -vf)^/^. When the photon reaches the far edge of the shell, the 
four- vector momentum becomes 



k' = ^E' 



/ 1 \ 

cosV^j 
sinV^j 




(17) 
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where ^2 is the angle between the normal vector of the void shell and the spatial three-vector of 
the momentum of the photon that leaves the void. 

As the photon leaves the shell, the four- vector momentum is converted to 

/ 1 \ 

smijj2 
\ J 

^ 7272[l + (^2-V2)COS'i/'2~V2V2] \ 
7272 [cos + (V2 - V2) - V2V2 COS ■lp2] 

simjj2 




~ R' ' 



\ 



) 



The velocities of the void are 



dr' 



R 



dt' 



(18) 



(19) 



(20) 
(21) 



where / = 1 , 2. From the connection conditions ([8]), and equation (flOl) . terms up to order 0((^^), 
V; can be calculated as 



1 



6(1 + ^0 



(22) 



where = ^(?;)and Shi = Snitd, = (l + SHi)^-(l+(l-^A(ti))5{tj)). Using the Friedmann equation, 
the parameters att = ti can be written in terms of those att = t2 as 



l + 5i H^fR',^^ 



i-^A(t2) 



1 + 62 Hf\R'J \\-nAmH2/Hif 



and 



K2 Hf \R\ 



(23) 



(24) 



From the geometry of the void in the internal comoving frame, at the order 0{^^), the relation 
between the void radius r,' j and r(,2 is approximately given by 



^vi sin «i<2 sin 1/^2 ■ 



(25) 
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The relation between the time t\ and ?2 can be written as 

"''^ dt' 



r' 1 cos ip'i + r[,2 cos ipj . 



(26) 



(27) 



The energy loss suffered between times t[ and t2 by a CMB photon that does not traverse the void 
is 

'EA ^ 1+Z2 

.-^1/ no void l+^l' 

where zi and Z2 are redshift parameters corresponding to ti and t2, respectively. The ratio of the 
temperature change for photons that traverse the void to that for photons that do not traverse the 
void is 

'E2\ l+zi 



AT 



1 / void 



1+Z2 



-1. 



(28) 



Equations (7-28) and can be solved recursively using ^2 as a small parameter. 
After a lengthy calculation, neglecting the terms of order 0(/t^ we find 

Ar 1 r q / 9 T. 

— = - COS iIj2(-26}j-5'h + (3 +45J6Hn^ 

+ (3 + 26,n)dH^m)COs2llj2 



(29) 



where 6h, S,„, (3, s, Q,„ are evaluated att = tx- For the matter-dominated EdS Universe, 6h = -Sm/3, 
17m = 1, and A = 0, we recover formula (34) from (Inoue & Silk 2006). Equation (|29l ) is the 
generalized formula for a spherical homogeneous dust-filled void in the flat-A FRW background. 
Plugging in P = for a compensating void (Sakai 1995; Sakai et al. 1999) B and the 

Hubble parameter contrast in linear order, 6h = ^,nSini^+ f~^(w))/2 (see the derivation and the 
definition of f{w) in the appendix), into equation (|29l ). the anisotropy formula for compensated 



^It turns out that the terms of order O(k^) are actually absent in AT /T if one writes 5h in terms of k and 5,„. The 
next leading order is O(k^). 

■'This approximation formula was originally derived for linear voids in low-density (flo < 1) FRW models. How- 
ever, from numerical analyses, it turns out that it can also be be applied to voids in flat-A FRW models if 6„, <C 1 
(Sakai 1995). 
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voids is 



1 



+ 



3 

54F2 



5Q 



2/3 



30ntl'F{ 1 + 2 cos(2^/;2)) + 9 (2^^"" + Bfi^/^ 



where 



5/6,1/3,11/6,-w 



(30) 



(31) 



and w is the effective equation-of-state parameter at time ?2- In contrast to the void model in the 
EdS background, there is an additional term which is linear in the matter density contrast 6,„. As 
shown in figure [H it can be interpreted as the contribution from the linear ISW effect, which comes 
from the change in the gravitational potential as the photon crosses the void. 

For a void at redshift z centered at an angular diameter distance d from the Earth, the angle 
tp2 between the light ray and the normal vector on the void surface at the time the photon leaves 
the shell is written in terms of the size/distance parameter b = R{z)ry/d and the angle 6 between 
the light ray and the line of sight to the void center as 

'ilj2{b,e)^^m\M/b). (32) 

Equations (l30l).(l3TI). and (|32|) give the anisotropy for compensating voids at arbitrary redshift z- 
For a void at z = in the flat-A FRW universe with Q„, o = 0.24, (Spergel et al. 2006) we have 
wq = ilmfl - 1 = -0.76. Then equations (l30l) and (|3TI) yield 

Ar 



0.087(cOS^^/^2)'^m,0 



- (0.021 +0.018 cos2V^2) 5; 



2 

m,0 



(33) 



where and 6^,0 are evaluated at t = to. Plugging ijj2 = (9 = 0) in equation (|30l) . we have the 
anisotropy 

AT 

— ~ (0.0875„,,o-0.0395^,o)Co (34) 
toward the center of the void. From equation (|34l) , we conclude that the contribution from the 
second-order effect, which is proportional to 5^ coherently enhances the contribution from the 



Here we do not call this the "Rees-Sciama effect" (Rees & Sciama 1968) because we deal with the relativistic 
second-order effect in which the curvature effect is no longer negligible. 
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linear ISW effect provided that the second-order correction for the void expansion and the wall 
expansion (which we discuss later) is negligible. In this case, we would observe a cold spot in the 
direction of the void. 

To generate an anisotropy Ar/T = 10"^, the void radius should be 

r,{z = 0) ^ 2 X 10'( 1 5^,0 1/0.2)-^/^ r^Mpc, (35) 

where H^^ = 3000 h~^Mpc. In comparison with the equivalent void model in the EdS background, 
the void radius can be smaller because of the coherent contribution owing to the linear ISW effect 
and the second-order effect. 

In order to study the second-order corrections, we introduce two nonlinear parameters (e, rj) 
that control the void and the wall expansion (see also Inoue & Silk 2006). These are defined as 

Oh = "mC>m - (^o„, (36) 

and 

P = -^Sm + v€- (37) 

If the matter inside the void near the boundary expands with the asymptotic velocity of the wall, 
i.e., 6h = P, then the non-linear parameter e can be written as 

e = ^ (nil' + 3fi^ - 66r]-5F-' ^ (38) 

Assuming 77 = 0, we have e > for fi^ > 0. For low-density universes = 0.2-0.3, we have 
e = (0.07-0.09)|(5|"^ As shown in figure[2l an increase in e leads to further redshifts of the photons 
due to reduction in the expansion rate inside the void. Therefore, for voids that expand with 
the asymptotic velocity, the second-order effect always enhances the linear ISW effect. On the 
contrary, for e < 0, the second-order effect leads to further blueshifts of the photons. Therefore, 
for voids that expand sufficiently faster than the asymptotic velocity of the wall (i.e. 6h > P), the 
second-order effect can reduce the redshift of photons due to the linear ISW effect. Furthermore, 
an increase in the velocity of the wall (i.e. 77 > 0) also leads to a suppression of the linear ISW effect 
because the photon is more Doppler blueshifted (see figure[2l right). Thus, the net redshift/blueshift 
of photons upon leaving the void depends on whether it is asymptotically evolving or not. In what 
follows, we define asymptotic voids as those with 6h = P and r] = 0. 



3. Origin of anomalies 



In order to explain the low-Z anomalies, we have proposed a model that consists of a pair of 
identical local voids in the direction (l,b) = (-30°, -30°) at z ~ 0.05 (Inoue & Silk 2006). In the 
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Fig. 1. — Amplitude of temperature anisotropy for a compensating void with (5„, = -0.3 as a func- 
tion of the mass density il,„ with the contribution from the term proportional to (5„, in equation (|30l) 
(solid curve), in equation (IB6I) (dashed curve) and that from the term proportional to 6f„ in equation 
(|30l) (dash-dotted curve). We assume ^2 = ''v2^2/^2 ' = 1 V'l = 0. 




Fig. 2. — Second-order effects. Temperature anisotropy profiles owing to the second-order effect 
are plotted as a function of the subtending angle 6. We assume a compensating void with a radius 
r,, = 200h~^ Mpc, Sm = -0.3, and the size/distance parameter b = 0.9. Le/?:(e,77)=(-0.3,0) (dotted 
curve), (0,0) (dashed curve), (0.3,0) (dash-dotted curve), and (0.6,0) (dash-dotted-dotted curve). 
Right:{e,ri)={0, 0.9) (dotted curve), (0,0.6) (dash-dot-dotted curve), (0,0.3) (dash-dotted curve), 
and (0,0) (dashed curve). The solid curve in the left panel denotes the temperature anisotropy 
profile owing to the linear ISW effect. 
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flat-A background, a similar model can work as well. For illustrative purpose, we have simulated 
skymaps for models that consist of identical asymptotic voids with a density contrast ~ -0.3 that 
are tangent to each other. Note that we are assumed to be just outside one of the voids. As the 
background cosmology, we have assumed a flat-A FRW model with i7,„o = 0.24. As shown in 
figure [3l the planar feature Z = 2, 3, 6 and the alignment feature / = 2, 3 can be well reproduced in 
our model that consists of identical voids with radius (2-3) x I0^h~^ Mpc. 

Although we still need an additional contribution that does not completely blur the fluctuation 
pattern by the voids, the chance of alignment between the quadrupole and octopole becomes large 
in our void model since otherwise no correlation between the two modes is expected. 

In order to study the large-angle correlation feature, we also calculated the two-point correla- 
tion function in real space C{6) (figure 4). For = 300/?"^ Mpc, the amplitude for the separation 
angle 9 > 30° is found to be C(9) ~ 100/iK^ which is marginally consistent with the observed 
WMAP 3-year values (see Spergel et al. 2006). 

Recently, it has been argued that the Shapley supercluster could be responsible for the anoma- 
lies on large angular scales (Cooray & Seto 2005; Rakic et al. 2006). However, as shown in figure 
[3l the amplitude in the direction to the Shapley supercluster (in the direction of the contact point of 
the two hypothesized voids) is somewhat smaller than those of the minima in either the quadrupole 
or the octopole. Furthermore, in the flat-A FRW model, the presence of the linear ISW effect 
(blueshift effect for mass concentration) can reduce the amplitude of the temperature fluctuations 
due to the Rees-Sciama effect (redshift effect for mass concentration) in the quasi-linear regime 
(locally Q,„ ~ 1). Because the highly non-linear structure (5 ^ 1) is concentrated in a particular 
region, it cannot strongly affect the large angular fluctuations. Thus, we conclude that such a sce- 
nario has a potential difficulty in explaining the anomalies in the FRW model with a cosmological 
constant. 



4. Mean contribution from dust-filled voids 

As shown in section 3, for linear voids \6,„ \ <^ 1 at z < 1, the temperature anisotropy due to a 
dust-filled void in the flat-A FRW background is proportional to (5,„rj^ due to the linear ISW effect. 
Assuming the cold dark matter (CDM) power spectrum, the mean amplitude of the temperature 
anisotropy due to a linear void is proportional to because |(5„,|rm.i oc for the linear regime. 
Therefore, we expect larger contributions from larger voids provided that |(5,„|r.„,..s ^ 1 (solid curve 
in the right panel in figure H]). For quasi-linear voids |5„,|r,«..s = O(10~'), or linear voids at high 
z, the second-order effect can also contribute to the anisotropy. Because we have |5„,|,..,„.i oc n}^^ 
for the non-linear regime, the contribution from the second-order effect is maximum at the scale 



-12- 




Fig. 3. — MoUweide projection maps for the cleaned CMB from the WMAP third- year data (Park 
et al. 2006) (left) and those for the temperature anisotropy owing to a pair of identical voids 
with radius = 300/i"^ Mpc (right) for modes 1 = 2-6. We assume a flat-A FRW universe with 
^m,o = 0.24. The boundaries of the void hemispheres are denoted by small light-gray disks. Note 
that the identical voids are tangential to each other and the apparently larger void is nearer to us. 
We also assume that we are outside the pair of voids. The size/distance ratios of the two voids 
are assumed to be & = 0.91 and b = 0.46. The separation angle of the void centers is 9s = 63" and 
the non-linear parameters are chosen to be (e,//) = (0.27,0), which correspond to an asymptotic 
void {5h = l3). The gray scale denotes the temperature fluctuations in which the maximum absolute 
value is set to 50^K. The north pole is aligned to the z-axis for which the angular dispersion of 
the quadrupole plus the octopole around the z-axis (de-Oliveira Costa et al. 2004) is maximal 
in the direction {l,b) = (-110°, 60°), which is close to the dipole direction at (/,Z?) = (-96°, 48°). 
The coordinate center is located at (/,fe) = (-30°, -30°). The Shapley supercluster is located in the 
direction of the point where the two voids contact each other (star). 



25 50 100 12b 1 = 175 



Fig. 4. — Two-point correlation function C(9) for the void model as in figure 3 (solid curve) and 
for the WMAP three year data (gray region) expected for a statistically isotropic sky with best-fit 
ACDM cosmology and for the best-fit flat ACDM model with = 0.04 and Q^fi = 0.24 (dash- 
dotted curve). The gray region denotes the maximum likelihood error at the 68% confidence level, 
which consists of the cosmic variance and the measurement errors. 
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r,, ~ 2 X I0^h~^ Mpc (dashed curve in the right panel in figure O. In the high redshift region 
z > 1, the contribution from the linear ISW effect is significantly reduced because the gravitational 
potential decay does not occur in the matter-dominant epoch. Similarly, the contribution from the 
second-order effect is also significantly reduced at z > 1 dH right). 

From these features, we conclude that local dust-filled quasi-linear voids with a radius r,, ~ 
200h~^ Mpc are natural candidates for explaining the origin of the claimed anomalies on large 
angular scales. 

However, one may argue that the possibility of having such large voids is extremely small. 
Indeed, if we assume a flat-A CDM cosmology, we find that a ~ 30% mass deficiency at the 
scale of = 200h~^ Mpc is only a 13cr result for linear fluctuations with ag = 0.74 or 1 la for the 
preWMAPS normalization = 0.9. B Assuming the CDM power spectrum (Bardeen et al. 1986), 
the expected temperature anisotropy for a linear (not asymptotic) void with radius = 200/7"' Mpc 
at z ~ is AT/T = 0(10"^). 

This inconsistency may be circumvented by considering a percolation process of small non- 
linear voids (Colombi et al. 2000; Hanami 2001). As A^-body simulations (Colberg et al. 2005) 
and local observations (Patiri et al. 2006) suggest, at low redshifts z < 1, non-linear voids with 
a typical radius of the order of O(10)/z"' Mpc cover almost the entire space. Therefore, small 
non-linear voids can percolate each other to make larger dust-filled local quasi-linear voids. In- 
terestingly, a low-density region with a 25% deficiency at z < 0.1 over a volume of 300/z"'Mpc^ 
has been discovered in the galaxy distribution in the south Galactic cap (Frith et al. 2003, 2005; 
Busswell et al. 2004). It seems that the underdense region consists of two local voids in the ranges 
0.03 < z < 0.06 and 0.07 < z < 0. 1, suggesting that the region consists of several dust-filled voids 
with a radius ~ 50h~^ Mpc. If such a low-density region can be regarded as a homogeneous 
dust-filled asymptotically evolving void, we expect an anisotropy in the local Hubble constant 
6h ~ -n^fSrn/e ~ 0.02. It can be even larger than 6h = {\ + f~\w))n,„S,n/2 ~ 0.04 if it is not 
asymptotically evolved. Although the substructure inside the void can also affect the anisotropy 
in the Hubble constant, we expect that the mean value over the low-density region will show a 
systematic deviation in comparison with the global value by several percent. Detection of a corre- 
lation between local large-scale structure and the anisotropy in the local Hubble constant would be 
strong evidence for the postulated local dust- filled large voids. 



^We consider a mean amplitude of linear fluctuations smoothed by a top-hat window function with a radius = 
200/j"' Mpc for cosmological parameters (f7„,,o, f^fc,o,/i, cr8,n) = (0.24,0.04,0.73,0.74, 1) 
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5. Summary and discussion 

In this paper, we have explored the large angular scale temperature anisotropy owing to the 
homogeneous local dust-filled voids in the flat-A FRW universe. A compensated asymptotically 
evolving dust-filled void with a comoving radius ry = (2-3) x lO^/i"^ Mpc and a density contrast 
6m = -0.3 will give a cold spot with anisotropy Ar/r~-10"^. A pair of such compensated local 
voids in the direction (/,&) = (-30°, -30°) might explain the quadrupole/octopole alignment and the 
planarity of quadrupole and octopole. As in the void model in the EdS background, the temperature 
anisotropy due to the dust-filled homogeneous void is proportional to (ry/H^^)^. For asymptotically 
evolving voids, contributions from the term that is proportional to 6^^ always enhance the linear 
ISW effect which is proportional to 6^- However, for local voids that expand sufficiently faster 
than the asymptotic velocity of the wall, contributions from the term oc 6f„ can suppress the linear 
ISW effect. Thus the enhancement/suppression depends on the non-linear dynamics of the void. 
We have also shown that the contribution to the CMB anisotropy from the linear ISW and the 
second order effects is more significant for local voids at z < 1 than those at high redshifts z > 1. 

Assuming the CDM power spectrum based on the standard inflationary scenario, the expected 
amplitude of the temperature anisotropy owing to the linear void with radius ry = 2x I0^h~^ Mpc 
is just AT/T = 0(10"^). However, at low redshifts z < 1, percolation of small non-linear voids 
could make large quasi-linear voids with a density contrast \5m\ ~ 0.1. The observed low-density 
region with a 25% deficiency at z < 0.1 over a volume of 300/z"^Mpc^ in the south Galactic cap 
(Busswell et al. 2004) might be evidence of such a structure. Detailed observational properties 
can be obtained by ongoing projects such as the 6dF galaxy survey (Heath et al. 2004). On 
the theoretical side, more elaborate modeling of the void percolation process is required because 
the signature on the CMB may depend on the configuration and the momentum of the substructure 
inside the void. Determination of non-linear parameters (e, 77) based on more realistic models is also 
an important task. Once these observational and theoretical issues are settled, we can determine 
whether new physics, such as primordial non-Gaussianity (Mathis et al. 2004) or an enhancement 
in the gravity in the dark energy sector (Farrar & Rosen 2006), may be required. 

Another implication unique to the dust-filled large voids postulated here concern increased 
dispersion in the locally measured Hubble constant as measured both in different directions and at 
different redshifts. As we have shown, for voids with a density contrast = -0.3, the expected 
fluctuation in the Hubble constant is as large as 2-4% (Tomita 2000, 2001). Although, in this 
paper, we have assumed that we are outside the dust-filled voids, it is also possible that we are 
inside the void. If the distance between the Milky Way and the wall of the void is sufficiently 
small, the ratio of our peculiar velocity to the light velocity will be v/c = 0(10"^) for a radius 
= (2-3) X 10^h~^ Mpc, and a density contrast 5,„ = -0.3. Therefore, the order of the predicted 
amplitude of the dipole is similar to that of the observed CMB dipole (Tomita 2003). In this case. 
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we expect a positively skewed dispersion in the Hubble constant 5h ~ 0.02-0.04 in the direction to 

the center of the void (approximately at il,b) = (-30°, -30°)). Accuracy measurement of ~ 1% in 
the Hubble constant would be achieved by increasing the number of samples of Type la supemovae 
from ~ 50 to 5000 provided that the error decreases as 1/ ^/N. 

Another approach to peculiar velocities uses galaxy surveys that apply different distance cali- 
brators (Pike & Hudson 2005; Watkins and Feldman 2007). The most optimistic individual galaxy 
distance estimates cited, using /-band surface brightness fluctuations, have an accuracy of 10% for 
half the sampled galaxies and extend to 40/2"' Mpc (Tonry et al. 2001). Other studies go deeper, 
for example the TuUy-Fisher relation has been applied to a sample of 3000 spiral galaxies out to 
200/i"'Mpc (Karachentsev et al. 2006) with distance errors of 15%. It should be possible to apply 
these techniques to test our hypothesis by targeting several thousand galaxies that sample the pe- 
culiar velocities induced by our nearest postulated voids and their associated shells, which extend 
over thousands of square degrees. 

Thus further study of correlation between the large scale structure and the variation in the 
locally measured Hubble constant is necessary in order to prove/disprove our void scenario. 

Our findings may shed new light on the issue of cross-correlation between the large-scale 
structure and the CMB. Within a flat ACDM model, recent works on the cross-correlation suggest 
a best-fit value Oa ~ 0.85 (Rassat et al. 2006; Cabre et al. 2006), which is larger than the best-fit 
value Qa = 0.70 from the joint analysis of the CMB and the large-scale structure (Tegmark et al. 
2004). Because the second order effect in quasi-linear asymptotic voids systematically enhances 
the linear ISW effect, the observed deviation could be explained by such quasi-linear or non-linear 
local voids in the flat-A universe with Qa = 0.70. 

If the low-£ CMB anomalies persist, our phenomenological model provides a simple astro- 
physical interpretation with no recourse to new physics. Perhaps the most promising tests of our 
hypothesis will come from improved sampling of Hubble flow anisotropics. If these directions 
coincide with the locations of our postulated voids, the case becomes more compelling. CMB 
polarisation maps and lensing studies could potentially reveal interesting signals. For example, the 
inflow pattern in the void wall will induce a small polarisation signal, as will the associated gravi- 
tational lensing of the CMB. The expected amplitude of the scattering angle of photons due to the 
postulated large local voids is ~ 5^^^ ~ 0.1°, corresponding to subdegree scales. These effects are 
small, amounting to an imprint on the ambient polarisation pattern of order a percent (~ v/c), but 
the phase structure would be unique and correlated with both the temperature map and the large- 
scale galaxy distribution. The global mean value of the Hubble constant will be slightly reduced, 
as the required dominance of voids biases the measured local Hq to be slightly high. Although 
the effect is small, it has potential significance for future dark energy surveys, which will combine 
CMB studies with deep galaxy redshift surveys. For example, the degeneracy between the Hubble 
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constant and the determination of the curvature of the universe using the CMB will need to be 
taken into account. 

One anomaly we have not sought to incorporate explicitly is the north-south power excess in 
the CMB. It may be tempting to attribute this to a directional deficiency of local power, as has been 
previously noted in deep galaxy counts (Busswell et al. 2004) and which in the present context 
consists of a compensated void complex that correlates the CMB and large-scale structure signals. 

The reduction in normalisation of erg to 0.74 makes it more difficult to account for large-scale 
features in the galaxy distribution without introducing extra large-scale power. Our example of 
tuned voids is one manifestation of extra power that would inevitably, if incorporated into the 
cosmological initial conditions of a simulation, boost the significance of large-scale features such 
as the SDSS "Great Wall" (Gott et al. 2005) and rare massive clusters (Baugh et al. 2004). 

We thank K. Tomita and N. Sakai for useful discussions and comments. This work is in part 
supported by a Grant-in- Aid for Young Scientists (17740159) from the Ministry of Education, 
Culture, Sports, Science and Technology in Japan. 



A. Linear perturbation of Hubble parameter 

In what follows, we calculate the Hubble parameter contrast at the matter- or A-dominant era 
in the linear order. Let us write the Friedmann equation at a point x in a perturbed universe in terms 
of the scale factor a, the Hubble parameter H, the matter density pm, the gravitational constant G, 
the curvature K, and the cosmological constant A as 

SttG ^ ^ K(x,t) A 

H\x,t) = -—p^(x,t) ^ + ^ (Al) 

3 a^{t) 3 

where the Hubble parameter H is defined on a comoving slice and the scale factor at the present 
time to is assumed to be ao = 1. In what follows, we introduce the total density p = Pm+ g^- The 
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linear perturbation component in equation (lAll) in the Fourier space for the flat FRW background 
can be written as 

2M//k=^Vk-^(^)Vk, (A2) 
where H is the spatially averaged Hubble parameter and the curvature perturbation IZ is defined as 

5K^ = ^k^TZk. (A3) 

If we can neglect the anisotropic stress, the relation between the Newtonian gravitational potential 
$k and the total energy density contrast 6^ is given by the Poisson equation 

kV 

- $k = 4^G(5pk. (A4) 
a J 

The time evolution of 6 is given by the continuity equation 

5p^ = -3(p+P)5Hk - SHSpi,, (A5) 

where p and = P are the spatially averaged total energy density and total pressure, respectively. 
The perturbed Raychaudhuri equation is (Liddle & Lyth 2000) 



6H^ = -2H6H^ ^ Vk + T - — (A6) 

3 3 \a/ p+P 

From equations (IA2HA6I) . we have 

^^-i<i'k+ ^^$k = -(1 +w)7^k, (A7) 

where the effective equation of state w = P/p is written in terms of the density parameter for the 
matter Q,„ o and that for the cosmological constant Qa,o at the present time to as 

w = ^- . (A8) 

^m.Q/(^A,0 a^)+l 



From the time derivative of the local Friedmann equation (|A2|) multiplied by a^, the continuity 
equation (IA5I ), and the Raychaudhuri equation (IA6I) . we have TZk = -H5P^/{p+P) = because the 
pressure perturbation vanishes SP^, = at the matter- or A- dominant era. Plugging $k = f{w)TZy^ 
into equation (IA7I) . since IZ^ = 0, /(w) satisfies 

2w(l+w)4^ + ^^^/+l + w = 0. (A9) 
dw 3 



The solution is 



fiw)=~{i+wy/\Fi 



5 1 11 

6'3'T'-" 



(AlO) 
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which satisfies the boundary condition / = -3/5 for w = 0. Plugging <l>k = f(w)7lk into equation 
(lA2l) and (lA4l) . we have 

5Hk \+f-'{w) 



H 

which yields 

l+f-\w) 



-4, (All) 



Sh = 2 ' (A 1 2) 

where the evolution of f(w) is given by equation (IA8l) and (lAlOl) . The equation (|A12|) is applicable 
to the matter- or A- dominated era at which the anisotropic stress is negligible. 



B. Linear integrated Sachs-Wolfe effect 

We evaluate the magnitude of the linear integrated Sachs-Wolfe (ISW) effect, which comes 
from the time evolution of the gravitational potential as the photon traverses the void. 

In terms of the gravitational potential \1/ , contribution from the ISW effect is approximately 
evaluated as 

(Bi) 

where At is the duration time for which the photons pass through the void. From the Poisson 
equation, in the Newtonian limit. At <^ H^^, the gravitational potential can be written in terms of 
the density contrast 6k = 6pk/ p inside the void in the Fourier space as 

= ~(2dd6k + d^Sk) (B2) 

where a dot means the partial time derivative • = ^ and k is the wavenumber. Plugging the peculiar 
velocity inside the void at comoving distance r from the center v = SnHar into the continuity 
equation 

— + a-W-\ = 0, (B3) 

ot 

we have 

d^5l = -3H^a^a5k, (B4) 

where a = (1 + f~^{w))/2. Writing the wavenumber in terms of the void comoving radius r,, as 
k ~ vr/r,,, equation (|B2I) and equation (IB4I) yield, 

Ar 1 



-e{2q + ^a)5, (B5) 

T Tl^ 
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where ^ = aVv/H'^ and q is the deceleration parameter. Thus, the anisotropy owing to the linear 
ISW effect is proportional to the density contrast 5 and In terms of the matter density contrast 
= Spm/ Pm ^nd the matter density parameter il,„, equation (IB5I ) can be written as 

— ~ S^^6m^,n, (B6) 

where 

1 /-„ 15 

s = — 



^ 3^]^----Q^;/V-M. (B7) 
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Fig. 5. — Mean amplitudes of temperature anisotropy for a linear void ^ 0.1 as a function 
of void radius at z = (/e/0,and as a function of void redshift z for = 200^"^ Mpc {right). 
Contribution from the linear ISW effect is denoted by solid curves and that from the second- 
order effect is denoted by dashed curves. We assume (e,0 = (0 0) and cosmological parameters 
{ftrnfi, ^b,o, h, as,n)^ (0.24, 0.04, 0.73, 0.74, 1) (Spergel et al. 2007). 



